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Abstract — Special attention has been paid to artificial 
groundwater recharge in water resource management in arid 
and semi-arid regions. Parameters considered in the selection of 
groundwater artificial recharge locations are diverse and 
complex. In this study, factors such as: slope, infiltration rate, 
depth to groundwater and electric conductivity (EC) are 
considered, to determine the areas most suitable for 
groundwater artificial recharge in aquifer in the Herat city of 
Afghanistan. Thematic layers for the above parameters were 
prepared, classified, weighted based on centroid method of 
de-fuzzification and integrated in a GIS environment by 
algebraic product operator of Fuzzy logic. Land use map of 
research area used to filter the artificial recharge map. The 
results of the study indicate that about 17.74% of the study area 
is suitable and 82.26% is unsuitable for artificial groundwater 
recharge. To validate the model a comparison between mean of 
water level points which located in suitable zone of suitability 
map and water level classes done. 

Index Terms — Artificial recharge, de-fuzzification method, 
GIS, Groundwater 


I. Introduction 

Water is indispensable for any life system to exist on earth 
and is a very important component for the development of any 
society [1]. A growing population and changing dietary trends 
mean a steeply rising water demand [2]. Demand for the 
worlds increasingly scarce water supply, is rising rapidly, 
challenging its availability for food production and putting 
global food security at risk, even as demand for water by all 
users grows, groundwater is being depleted [3]. Groundwater 
as a source of water supply has great advantages over surface 
water from streams, rivers, or lakes [4] and with increasing 
demands for water, ground resources are gaining much 
attention [5]. 

Afghanistan’s climate is arid to semi-arid where the 
weather is cold in winter and hot and dry in summer with 
temperature that ranges from -20°C in winter to 50°C in 
summer. The annual precipitation in Afghanistan varies 
according to region, ranging from75mm in the southwest up to 
1270 mm in the northeast with an average annual precipitation 
of 300mm [6]. According to the Ministry of Energy and Water 
[7] of Afghanistan, long-term total annual precipitation in 
Afghanistan is 164 billion m 3 , evaporation is 87 billion m 3 
(53% of the total annual precipitation), surface water runoff is 
61 billion m 3 (37% of the total annual precipitation) 
groundwater recharge is 16 billion m 3 (10 % of the total 
annual precipitation), and the total annual available water 
resources are 77 billion m 3 (Table 1). 

Nasir Ahmad Gesim, Computer Science & Intelligent systems, 
University of the Ryukyus, Okinawa, Japan 

Takeo Okazaki, Computer Science & Intelligent systems, University of 
the Ryukyus, Okinawa, Japan, +81-98-895-8903 


II. Research propose 

Above statistics indicate that just 10 % of the total annual 
precipitation recharge groundwater naturally and since more 
than 15% of the agricultural land is being irrigated by 
groundwater in karez wells fed by springs and shallow wells 
in Afghanistan [7] therefore, it is important to augment the 
groundwater resource by artificial recharge. Effective 
management of aquifer recharge is becoming an increasingly 
important aspect of water resource management strategies [8]. 
Conservation of soil and its proper utilization must also be 
considered as a natural resource, in water resource 
management plans. So, artificial recharge of groundwater will 
be helpful to rising water level. 

III. RELATIONAL STUDIES 

There are many approaches for selection of suitable 
location for artificial recharge (AR). The application of 
traditional data processing methods in site selection for 
artificial groundwater recharge is very difficult and time 
consuming, because the data is massive and usually needs to 
be integrated [11]. Different operators of Fuzzy logic such as 
AND, OR, algebraic sum, algebraic product and gamma are 
capable to develop information in different thematic layers 
and integrating them with sufficient accuracy and within a 
short period of time. The application of these methods is 
indispensable for such analyses. Many studies have used 
fuzzy operators for locating most suitable sites for artificial 


Table .1 Annual precipitation data by basin in 
Afghanistan (MEW, 2016) 


Basin 

Area (km 2 ) 

Annual mean 
precipitation 
(mm) 

Total 

precipitation 
(billion m 3 /year) 

% of total 
precipitation 

Kabul 

108392 

298 

32.3 

20 

Helmand 

202,006 

180 

36.3 

22 

Amu 

101,498 

393 

39.8 

24 

North 

78,099 

268 

20.9 

13 

Harirod 

162,659 

210 

34.1 

21 

Total 

652,654 

270 

163.3 

100 


recharge. (Nirmala et al,2011) [9] used hybrid algorithm to 
stud artificial recharge of groundwater in Sathyamangalam 
and Melur villages, Chennai. (Ghayomian et al, 2007) [11] 
applied Fuzzy logic among GIS techniques to determine most 
suitable areas for artificial groundwater recharge in a 
coastalaquifer in Gavbandi Drainage Basin, Also, (Nouri et 
al,2005) [12] have had the same studies and used Fuzzy 
algebraic product to carry out their study. 
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IV. Research area 

The research area (Herat city) located in the center of Herat 
province, Afghanistan located between (34.248° and 34.474°) 
latitudes, and (61.942° and 62.442°) longitudes with an area 
of 730 km 2 (Fig. 1). Land-surface elevation in the research 
area ranges from 858 to 1636 above sea level, with an average 
of 1247m. The mean annual precipitation is recorded as 
210mm [7]. Exploitation of groundwater resources in the 
study area includes use of qanat, springs, and deep and 
semi-deep wells. The average well discharge is approximately 
200 Lit/min. The research area also consists of 230 wells 
where water is withdrawn from the alluvial fan and the well 
depths range between 7 and 90m. The general trend of 
groundwater flow is from the east to west. 

V. Materials and methods 

A. Thematic layers 

Among effective factors in locating suitable areas for 
artificial recharge slop, depth to groundwater, infiltration rate, 
electric conductivity and land-use factors were selected and 
examined [13]. Among these, slop, depth to groundwater, 
infiltration rate, electric conductivity factors which have 
direct impact on artificial recharge were used as principle 
factors and land use factor which shows feasibility of 
implementation were used as a filter. 

The main source of groundwater is the water derived from 
rain and snow-melt that has permeated through the alluvium 
and usually this water has a good quality. After permeating the 
soil, rain water quality will change because of the contact with 
alluvium and dissolving different minerals in it. This 
qualitative change depends on constituting particles of 
aquifer, duration of contact with the bedrock, utilization rate 
of the groundwater, and groundwater level. Therefore, the 
quality of groundwater in alluvium as an essential parameter 
was investigated in the artificial recharge. Since electric 
conductivity and total dissolved solids (TDS) variations have 
similar trends, the EC factor was used as an indicator for 
water quality. Fig.3a illustrate electric conductivity map of 
research area which prepared based on the actual data of 169 
well and Inverse Distance Weighting (IDW) method 
interpolation [24]. 

Infiltration in its most narrow and precise sense can be 
defined as the process water entering into soil through the soil 
surface. Although a distinction is made between infiltration 
and percolation the movement of water within the soil the two 
phenomena are closely related since infiltration cannot 
continue unimpeded unless percolation removes infiltrated 
water from the surface soil. The soil is permeated by 
noncapillary channel through which gravity water flows 
downward towards the ground water, following the path of 
least resistance. Capillary forces continuously divert gravity 
water into pore spaces, so that the quantity of gravity water 
passing successively lower horizons is steadily diminished. 
This leads to increasing resistance to gravity flow in the 
surface layer and a decreasing rate of infiltration as a storm 
progresses. The rate of infiltration in the early phases of a 
storm is less if the capillary pores are filed from a previous 
storm. There is maximum rate at which water can enter soil at 
a point under a given set of conditions; this rate is called the 
infiltration capacity. The actual infiltration rate equals the 


infiltration capacity only when the supply rate rainfall 
intensity less rate of retention) equals or exceeds. Fig.3b 


Herat Province 

Afghanistan 



f ■* 

■ 

t JP 

/fn 



e2 c 0'0"E 62°10'0“E 

62°20'0'E 62°30'0“E 



62WE 62°10'0'E 62°20'0'E 

Fig. 1 Research area 

62°30’0-E 

Table .2 Membership functions of thematic layers 

Thematic layers 

Classes 

membership 


0-2 

0.83 


2-4 

0.69 

slope(%) 

4-8 

0.44 


8-10 

0.18 


>10 

0.01 


13 

0.42 

Infiltration Rate(mm/hr) 

18 

0.49 

25.9 

0.58 


61.2 

0.78 


0-4 

0.01 


4-15 

0.42 

Water Level(m) 

16-25 

0.65 

26-35 

0.7 


36-45 

0.81 


46-63 

0.88 


740-1300 

0.89 


1400-1900 

0.82 

Electric Conductivity 
(ms/s) 

2000-2400 

0.74 

2500-3000 

0.65 


>3000 

0.27 
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illustrates the infiltration rate of research area based on the 
Herat city soil types and hydrologic soil properties [19]. 

The water depth in infiltration basins should select 
carefully, while high hydraulic heads produced by deep water 
result in high infiltration rates, they also tend to compress 
clogging layers. Fig.3c shows the water level situation of 
Herat city based on 230 actual data that collected during 
overseas research. 

Another extremely important factor for identification of 
suitable site is slope. This parameter plays an important role 
to control factors like runoff, erosion, material transportation 
and permeability. Fig. 3d shows slope percent in research area 
which classified in 5 classes. 

B. Fuzzy logic 

In the classical setting an element either belongs to a set 
or not. If A is a classical set, then the formula X G A is either 
absolutely true or absolutely false. In the case of a fuzzy set A, 
an element x can attain more than two degrees of its 
membership. Thus, the formula X G A may be only partially 
satisfied [13]. Hence then, fuzzy subsets have been applied to 
diverse field [14]. In contrast to Boolean logic; no certainty 
exists in fuzzy logic. Therefore, no unit is satisfactory or 
unsatisfactory in this logic. Since in previous studies on 
artificial recharge of groundwater membership values of 
different classes of a thematic layer have assigned empirically 
so, here tried to define these values using fuzzy logic. 

Given two or more maps with fuzzy membership functions 
for the same set, a variety of operations can be employed to 
combine the membership values together [11]. Zimmermann 
and Zysno (1980) discuss a variety of combination rules [15]. 
There are five operators that were found to be useful for 
combining exploration datasets, namely the fuzzy AND, 
fuzzy OR, fuzzy algebraic product, fuzzy algebraic sum and 
fuzzy gamma operator [10]. The fuzzy OR and AND 
operators are used, only one of the contributing fuzzy sets 
influences the resultant value. The fuzzy algebraic sum and 
algebraic product operators make the resultant set larger than, 
or equal to the maximum value and smaller than, or equal to 
the minimum value among all fuzzy sets, respectively. The 
fuzzy gamma operator has the value between that of the fuzzy 
algebraic product operator and that of the fuzzy algebraic sum 
operator. Fuzzy algebraic product which selected here 
because of its high sensitivity in specifying artificial recharge 
areas [11] defined as [14]: 

Mcomb = H/', (1) 

where jil i is the fuzzy membership function for the ith map, 

and i = (1, 2, 3,...,ft), maps are to be combined. In fuzzy 
algebraic product operator as a t- norm, the weight of 
compositional layer in the multilayer intersection section is 
equal to their products and for other sections is zero. 
Therefore, mentioned operator has a decrease effect [15]. 
Although the fuzzy algebraic product gives an output that is 
decreasive in nature, it does utilize every membership value to 
produce the result, unlike the fuzzy minimum [16]. Fuzzy 
logic with range of zero to one is considered for different 
satisfactory levels. 


In this study as a new approach centroid method 
de-fuzzification which showed better result than other method 
of de-fuzzification [18] used to determine membership degree 
of each class of thematic layers. Basically, this logic system 
consists of the following: 

1. Fuzzification: Converting the crisp inputs to 

membership functions which comply with intuitive 
perception of system status. 

2. Rules Processing: Calculating the response from system 

status inputs according to the pre-defined rules matrix 
(control algorithm implementation). 

3. Inference: Evaluating each case for all fuzzy rules 

4. Composition: Combining information from rules 

5. De-Fuzzification: Converting the result to crisp values. 
To prepare standard fuzzy set as an input Gaussian Function 
used to convert crisp values 

(■ X-b ) 2 

fix) = ae 2 ° 2 (2) 

Where the parameter a is the height of the curve's peak (here 
in fuzzy logic is 1), b is the position of the center of the peak 
and c the standard deviation . After processing input in the 
inference engine as illustrated in Fig.2, the result of the 
process must be converted to a crisp value, by composition of 
the crisp values the membership curves for different classes of 
thematic layers will prepared. 

VI. Experiments 

The objective of this study is to firstly determine 
membership values for different classes of four mentioned 
layers and secondly to integrate thematic layers in order to 
prepare groundwater artificial recharge map of research area. 

Table.2 shows membership degrees for different classes of 
thematic layers which prepared by centroid method of 
de-fuzzification in Fuzzy Inference System (Fig.2). Then 
algebraic product operator of fuzzy logic was used to overlay 
these layers. 

Suitability map of artificial recharge which developed by 
applying fuzzy logic models to the thematic layers shows that 
20.08% percent of area is suitable (Fig.4e). Land-use map of 
research area as illustrated in figure.4f classified in five 
classes and coded as one (suitable) and zero (unsuitable). This 
classification is applied to the map of areas suitable for 
recharge, as a filter. After filtration the suitable area decreased 
to 17.74% as illustrated in Fig .4g 


Fuzzifier ^EnTne^ Defuzzifier 
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Fig. 3 thematic layers (a) electric Conductivity, (b) 
infiltration rate, (c) depth to groundwater, (d) slope 


To validate the model, a comparison between mean of water 
level points which located in suitable zone of Fig.4g and water 
level classes of Fig.3c done [22]. The average of water level 
points which located in suitable zone is 10.84m and since non- 
accepted range of water level is just first class of Fig.3c where 
depth to groundwater is less than 4m as its membership degree 
illustrated 0.01 in Table.2. Thus, it can be inferred that the 
artificial recharge zones delineated by remote sensing and 
GIS techniques are reliable. 

VII. Conclusion 

Groundwater resources provide the most important amount 
of water resources, and many agricultural units have been 
subjected to negative balance of water due to overharvesting 
of wells in Afghanistan [7]. Hence proper management of 
groundwater resources is necessary and artificial recharge of 
groundwater can be helpful in management of groundwater 
and increasing the groundwater level. One of the most 
important factors in successful recharge of groundwater 
resources is locating suitable areas for these projects. 
Therefore, site selection of artificial recharge suitable areas is 
very important. 

Four factors namely, depth to groundwater, electric 
conductivity, slope, and infiltration rate parameters were 
explored, classified, weighted and overlaid. To overlay these 
layers, fuzzy logic was used, and the study areas were divided 
into two suitable and unsuitable classes for artificial recharge. 
Results showed that based on overlay fuzzy logic 17.74% of 
study area was suitable for artificial recharge. In this matter, it 
can be said that suitable areas have been decreased from 
20.08% to 17.74% using land use filtering. In fact, land use is 
one of the major factors of water resources management 
restriction. This artificial recharge zone map can serve as a 
guideline for the planners/decision-makers as well as 
practicing hydrogeologists for the sustainable management of 
groundwater resources in the study area. 
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Fig.4 (e) Suitability map of AR before filtration, (f) land 
use map, (g) suitability map of AR after 
Filtration 
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